Generation of quasi-periodic waves and flows in the solar 
atmosphere by oscillatory reconnection 
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ABSTRACT 

We investigate the long-term evolution of an initially buoyant magnetic flux tube emerging into 
a gravitationally-stratified coronal hole environment and report on the resulting oscillations and 
outflows. We perform 2.5D nonlinear numerical simulations, generalizing the models of McLaugh- 
lin et al. ( 2009) and Murray et al. ( 2009). We find that the physical mechanism of oscillatory 
reconnection naturally generates quasi-periodic vertical outflows, with a transverse/swaying as- 
pect. The vertical outflows consist of both a periodic aspect and evidence of a positively-directed 
flow. The speed of the vertical outflow (20 — 60 km/s) is comparable to those reported in the 
observational literature. We also perform a parametric study varying the magnetic strength of 
the buoyant flux tube and find a range of associated periodicities: 1.75 — 3.5 min. Thus, the 
mechanism of oscillatory reconnection may provide a physical explanation to some of the high- 
speed, quasi-periodic, transverse outflows/jets recently reported by a multitude of authors and 
instruments. 



Subject headings: Magnetic Reconnection — MHD 
— Sun: oscillations 

1. Introduction 



Improvements in the spatial and temporal res- 
olution of solar observations have led to a recent 
deluge in reported magnetohydrodynamic (MHD) 
wave motions (e.g. see Nakariakov & Verwichte 
2005; De Moortel 2005; Banerjee et al. 2007; 
Ruderman & Erdelyi 2009; Goossens et al. 2011; 
McLaughlin et al. 2011 for a recent list). Here, 
we focus on observations of transverse motions in 
the solar atmosphere (e.g. Tomczyk et al. 2007 
De Pontieu et al. 2007; Cirtain et al. 2007 
Erdelyi & Taroyan 2008; He et al. 2009a; 2009b 
Liu et al. 2009; 2011; Morton et al. 2011 
Okamoto et al. 2011). These transverse motions 



Waves — Sun: activity — Sun: magnetic fields 

have been called Alfven waves by some authors, 
although this is subject to debate and they are al- 
ternatively interpreted as kink waves (e.g. see ar- 
guments by Erdelyi & Fedun 2007; Van Doorsse- 
laere et al. 2008). The dispute rests not with the 
observations themselves, but with the appropri- 
ate interpretation: MHD wave modes of an over- 
dense cylinder versus MHD waves of a homoge- 
neous plasma. These arguments and others, e.g. 
whether or not a stable waveguide actually exists 
in the solar atmosphere, are not the focus of this 
current paper. 

Tomczyk et al. ( 2007) utilised the CoMP /Coronal 
Multi- channel Polarimeter instrument to report 



on ubiquitous, small-amplitude, transverse dis- 
turbances, propagating along magnetic field lines. 
The authors do not report on how the waves are 
generated, but do speculate the waves may origi- 
nate from within the chromospheric network that 
forms the footpoints of the observed loops. Dc 
Pontieu et al. ( 2007) used Hinode/SOT measure- 
ments in an attempt to reveal Alfven/transversal 
waves in the chromosphere with strong amplitudes 
(10 — 30 km/s) and periods 100 — 500 seconds. Ca 
II H-line images also reveal a plethora of dynamic, 
jet-like extrusions called chromosphere spicules, or 
type II spicules. These spicules undergo a sway- 
ing/oscillatory motion perpendicular to their own 
axis, which the authors described as Alfvenic mo- 
tions. Again, this interpretation is disputed by 
other authors (e.g. He et al. 2009b; Verth et al. 
2011) who interpreted these spicule oscillations as 
kink waves, due to the fact that spicules are over- 
dense in comparison with the ambient plasma. 

De Pontieu et al. ( 2011) report on a link be- 
tween chromospheric spicules/jets and their coro- 
nal spicules/jets counterparts, i.e. suggesting a 
mechanism for imparting chromospheric plasma 
into the corona. The coronal spicules are strongly 
heated and are seen to rapidly propagate upwards, 
but the authors report that there are currently no 
models for what drives and heats the observed jets 
(see also a review by Sterling 2000). Okamoto & 
De Pontieu ( 2011) report on the statistical prop- 
erties of transverse (Alfvenic) waves along spicules 
and report median velocity amplitudes and peri- 
ods of 7.4 km/s and 45 seconds, respectively (see 
also review by Zaqarashvili & Erdelyi 2009). 
Mcintosh et al. (2011) reported transition re- 
gion observations of ubiquitous, transverse (sway- 
ing/Alfvenic) motions that are outwardly propa- 
gating, with amplitudes ~ 20 km/s and periods 
<~ 100 — 500 seconds, energetic enough to heat the 
fast solar wind. Again, the authors note that the 
challenge remains to understand how and where 
these waves are generated in the solar atmosphere. 

Thus, transverse/swaying motions have been 
observed over a range of temperatures, wave- 
lengths, speeds and scales. However, the origin of 
these propagating, transverse oscillations remains 
a mystery. Liu et al. ( 2011) summarise possible 
generation mechanisms including an oscillating 
wake from a CME or periodic reconnection (e.g. 
Chen & Priest 2006; Sych et al. 2009). 



1.1. Waves versus Flows Interpretation 

In addition to transverse MHD wave observa- 
tions, slow MHD waves have also been observed 
in the solar atmosphere (e.g. Ofman et al. 1997; 
DcForest & Gurman 1998; Berghmans & Clette 
1999; De Moortcl et al. 2000; Ofman & Wang 
2008; Erdelyi & Taroyan 2008). More recently, 
this propagating slow wave interpretation has been 
challenged. De Pontieu & Mcintosh ( 2010) show 
that Hinode/EIS measurements of intensity and 
velocity oscillations of coronal lines are driven by 
a quasi-periodically varying component in the blue 
wing of the emission line, i.e. providing evidence 
of quasi-periodic upflows. Moreover, such upflows 
(<~ 50 — 150 km/s in the line-of-sight) have also 
been reported in coronal loop footpoints in active 
regions (De Pontieu et al. 2009) and quiet-Sun 
regions (Mcintosh & De Pontieu 2009a). These 
authors also suggest a direct link between these 
outflows/propagating disturbances and chromo- 
spheric spicules, i.e. the fountain-like jets that 
protrude into the corona (Mcintosh & De Pontieu 
2009b). 

Mcintosh et al. ( 2010) analysed STEREO 
observations of polar plumes and identified high- 
speed upflows, again as opposed to previous in- 
terpretations as propagating slow waves. These 
authors observed high-speed jets travelling along 
the plumes with mean velocity of 135 km/s and re- 
peating quasi-periodically (repeat times 5-25 min- 
utes). Using SDO/AIA to extend this analysis to 
plume-like structures originating from equatorial 
coronal holes and quiet-Sun regions, Tian et al. ( 
2011) found that the outflows are not restricted to 
plumes. These authors reported that the outflows 
exhibit transverse, swaying motions and probably 
originate from the magnetic network of the quiet- 
Sun and coronal holes. In contrast, Verwichte et 
al. ( 2010) argue that a slow mode interpreta- 
tion remains a valid explanation for the observed 
quasi-periodic intensity perturbations, and show 
that slow waves inherently have a bias towards 
emission in the blue wing of the emission line due 
to the in-phase nature of the velocity and density 
perturbations. 

Thus, it is clear that there is a need for a bet- 
ter understanding of the physical generation mech- 
anism for these ubiquitous, transverse, propagat- 
ing motions, and that both quasi-periodic motions 



and possible flows must be kept in mind. This pa- 
per attempts to describe such a physical mecha- 
nism by utilising a dynamic, reconnection numer- 
ical model. 

1.2. Flux Emergence &: Oscillatory Recon- 
nection 

Magnetic field is continuously emerging on the 
Sun over a range of scales (see Archontis 2008 
for a comprehensive review). Magnetic flux tubes, 
formed at the tachocline, rise buoyantly through 
the convection zone. As they reach the photo- 
sphere, their buoyant-rise ends and an instability 
allows the magnetic flux to penetrate into the so- 
lar atmosphere. The subsequent evolution of the 
newly-emerged flux is then dominated both by the 
properties of the rising flux tube itself and the 
pre-existing magnetic topology the tube emerges 
into. Flux emergence is well described in the ex- 
isting literature, as is the dynamic reconnection 
associated with the collision of newly-emerged flux 
and pre-existing magnetic field (e.g. Shibata et al. 
1992; Yokayoma & Shibata 1995; 1996; Archon- 
tis et al. 2004; 2005; 2006; 2007; Isobe et al. 
2005; 2006; Murray et al. 2006; Galsgaard et al. 
2007; Moreno-Inserti et al. 2008, and references 
therein). 

Shibata et al. ( 2007) reported Hinode/SOT 
observations of the ubiquitous presence of chro- 
mosphcric anemone jets (velocity 10 — 20 km/s). 
These numerous, small-scale jets, seen in, e.g., Ca 
II H broadband, display an inverted Y-shape, i.e. 
the characteristic shape of anemone jets (e.g. Shi- 
bata et al. 1994; Yokoyama & Shibata 1995). 
The anemone shape is formed as a result of mag- 
netic reconnection between an emerging magnetic 
bipole and a pre-existing vertical field. 

Reconnection can occur when strong currents 
cause the magnetic ficldlines to diffuse through 
the plasma and change the connectivity (Parker 
1957; Sweet 1958; Petsheck 1964). In 2D, re- 
connection can only occur at null points (Priest 
& Forbes 2000). Dungey ( 1953) reported that a 
perturbed X-point can collapse if the footpoints 
are free to move, Mellor et al. ( 2002) stud- 
ied the linear collapse of a 2D null point, and 
Imshennik & Syrovatsky ( 1967) described the col- 
lapse with an exact, nonlinear solution of the ideal 
MHD equations. However, these papers do not in- 
clude the effect of gas pressure, which acts to limit 



the growth of the current density. In considering 
the relaxation of a 2D X-typc null point, Craig 
& McClymont ( 1991) found that free magnetic 
energy is dissipated by the phenomenon of oscil- 
latory reconnection , which couples resistive diffu- 
sion at the null to global advection of the outer 
field. McLaughlin et al. ( 2009) investigated the 
behavior of nonlinear fast magnetoacoustic waves 
near a 2D X-point and found that the incoming 
wave deforms the null point into cusp-like point 
which in turn collapses to a current sheet. The 
system then evolves periodically through a series 
of horizontal/vertical current sheets with associ- 
ated changes in connectivity, i.e. the system dis- 
plays oscillatory reconnection. Longcope & Priest 
( 2007) investigated the diffusion at the null of a 
2D current sheet subjected to a suddenly enhanced 
resistivity, finding that the diffusion couples to a 
fast mode which propagates the current away at 
the local Alfven speed. 

Of particular importance to the work presented 
in this paper is that of Murray et al. ( 2009). Flux 
emerging into a pre-existing field has been studied 
in great detail before, but Murray et al. ( 2009) 
were the first to investigate the long-term evolu- 
tion of such a system, i.e. many previous sim- 
ulations end once reconnection is first initiated. 
Murray et al. utilised a stratified atmosphere per- 
meated by a unipolar magnetic field (representing 
a coronal hole) and investigated the emergence of 
a buoyant flux tube. Murray et al. found that 
a series of 'reconnection reversals' take place as 
the system searches for equilibrium, i.e. a cycle of 
inflow/outflow bursts followed by outflow/inflow 
bursts. Thus, the system demonstrates oscillatory 
reconnection (e.g. Craig & McClymont 1991; 
McLaughlin et al. 2009), initiated in a self- 
consistent manner. Murray et al. also detail the 
physics behind the phenomena. 

The aim of this paper is to further generalize 
the model of Murray et al. ( 2009) and to detail 
the oscillatory outputs due to oscillatory recon- 
nection (Craig & McClymont 1991; McLaughlin 
et al. 2009). We will also investigate the depen- 
dency and robustness of the model by varying the 
initial magnetic strength of the buoyant flux tube. 

The paper has the following outline: the numer- 
ical model is detailed in §2, brief recall of Murray 
et al. ( 2009) is described in §3 and the quasi- 
periodic outputs are reported in §3.1. §4 investi- 
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Fig. 1. — (a) Equilibrium magnetic field: vertical magnetic field representing a unipolar coronal hole and 
curved fieldlines of our initial flux tube, centered at (x,y) = (0, —4.25) Mm. (b) Equilibrium conditions of 
the numerical model: temperature (solid), gas pressure (dashed), density (dot-dashed). In both figures, the 
solar interior (SI), photosphere (PH), transition region (TR) and corona (CR) are indicated by the dotted 
lines. Only a subset of the full numerical domain is shown. 



gates the dependency of the model to the initial 
strength of the buoyant flux tube and conclusions 
arc presented in §5. 

2. Numerical Model 

We consider the two-dimensional, nonlinear, 
compressible, resistive MHD equations, including 
gravitational effects: 
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where p is the mass density, v is the plasma ve- 
locity, B the magnetic induction (usually called 
the magnetic field), p is the plasma pressure, 
p = 4-7T x 10 _7 Hm _1 is the magnetic permeabil- 
ity, acceleration due to gravity g = —gy, a is the 
electrical conductivity, p = 1/ pa is the magnetic 
diffusivity, e = p/p (7 — 1) is the specific inter- 
nal energy density, 7 = 5/3 is the ratio of specific 
heats and j = V x B//x is the electric current den- 
sity. 

We solve these governing equations using the 
LARE2D numerical code (Arber et al. 2001) 
which utilizes artificial shock viscosity to introduce 
dissipation at steep gradients, and the details of 



this technique, often called Wilkins viscosity, can 
be found in Wilkins ( 1980). Thus, Q s hock rep- 
resents the viscous heating at shocks. Heat con- 
duction and radiative effects are neglected in the 
present study. 

We now introduce a change of scale to non- 
dimensionalise all variables. Letting v = vov*, 
B = SB*, x = Lx* , y = Ly* , p — pop*, 
p = pop*, V = V*/X, t = t Q t*, A = BLA*, 
9 = .9o and r/ = rjo, where * denotes a dimen- 
sionlcss quantity and vo, -B, L, po, po, to, go 
and po are constants with the dimensions of the 
variable they are scaling. Here, A = A z z is 
the vector potential. We then set B/^/ppo = vo 
and vo = L/to (this sets vo as a constant back- 
ground Alfvcn speed). We also set ijoto/L 2 = -R" 1 , 
where R m is the magnetic Reynolds number. This 
process non-dimensionalises equations (1) and un- 
der these scalings, t* = 1 (for example) refers to 
t = to = L/vo; i.e. the time taken to travel a 
distance L at the background Alfven speed. For 
the rest of this paper, we drop the star indices; the 
fact that all variables are now non-dimensionalised 
is understood. 

The values returned from equations (1) arc 
made dimensional using the following choice of so- 
lar constants: (photospheric) pressure scale height 
L = 1.7 x 10 5 m, time to = 25 seconds, veloc- 
ity v = L/to = 6.8 x 10 3 m/s, density p = 
3.0 x 10~ 4 kg/m 3 , pressure po = 1.2 x 10 4 Pa, 
temperature To = 5.6 x 10 3 K, magnetic field 
B = 1.3 X 10 3 G and gravity 50 = 270 m/s. 
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Fig. 2. — Contours of current density (|J|, units Am -2 ) and selection of fieldlines at times t = 30, 33, 36, 
39, 42 and 45mins. Dotted lines denote horizontal layers of the equilibrium solar atmosphere. Note that the 
numerical domain plotted and colour bars change between the two rows. 



2.1. Initial and Boundary Conditions 

To simulate flux emerging into the solar atmo- 
sphere, we follow the numerical set-up of Murray 
et al. ( 2009). We consider a numerical domain 
comprising of four horizontal layers (Figure la). 
Above the lower boundary, we have a solar interior 
that is marginally stable to convection, a 5600 K 
isothermal photosphere, a transition region with a 
rapid (power law) temperature increase and, above 
this, an isothermal corona. Each of the layers 
is initially in hydrostatic equilibrium. Figure lb 
shows the equilibrium conditions. 

To model a unipolar corona hole, we chose an 
initial magnetic field to be a negatively-directed 
vertical magnetic field of 19.5 G. The coronal hole 
temperature, density and field strength are taken 
from Zhang et al. ( 2007) and Baker et al. ( 2009). 

To initiate flux emergence, a magnetic flux tube 
is placed in the solar interior at x — and a 
depth of y = —1.7 Mm. In cylindrical coordi- 
nates, the magnetic flux tube is chosen to have 
B = Bo (o,are- r2 / R2 ,e~ r2 / R2 y We choose B = 



3.25xl0 3 G, R = 4.25xl0 5 m, and a = -0.064x2tt 
for each L length in the axial direction. The 
buried flux tube is set in radial force balance and 
in thermal equilibrium with the external plasma. 
Thus, a density difference exists such that the flux 
tube is buoyant relative to its surroundings and, 
at the start of the simulation, will rise bodily up- 
wards/towards the photosphere. 

We chose a numerical domain — 68 < x < 
68 Mm, -4.25 < y < 63.75 Mm using 1600 x 832 
gridpoints of uniform spacing. All boundaries are 
fixed in x and y. Convergence testing was carried 
out with double and half the resolution. 

3. Flux Emergence &: Recreation of Mur- 
ray et al. results 

Flux emergence is well documented in existing 
literature (e.g. Murray et al. 2006, Archontis 
et al. 2008, as well as the papers listed in §1.2 
and references therein). Figure 2 (top row) shows 
that the buoyant magnetic tube rises and emerges 
into the model atmosphere, and that the emerging 
flux compresses the pre-existing magnetic field as 
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Fig. 3. — Contours of current density (|J|, units A m -2 ) and selection of fieldlines at times t = 46.7, 47.8, 
49.2 and 50.8 mins. Dotted lines denote a horizontal layer of our equilibrium solar atmosphere (i.e. change 
from transition region to coronal temperature profile). Note that we plot different axes compared to Figure 
2. 



it expands. To the north-west side of the emerging 
flux, the magnetic field is directed positively out of 
the solar surface whereas the neighboring coronal 
hole is directed in the opposite direction. Thus, a 
current sheet builds up at this interface between 
the two flux systems and this can be clearly seen 
in the bottom row of Figure 2. Reconnection com- 
mences at t = 41.6 mins. 

The system demonstrates the phenomenon of 
oscillatory reconnection as it searches for an equi- 
librium. In Figure 3, we see that at t = 46.7 mins 
the system forms a current sheet located around 
-3.5 < x < -2.5 Mm, 4.5 < y < -5.5 Mm, at 
an angle of approximately 7r/4 relative to the pos- 
itive x— direction. In this paper, we shall refer to 
a current sheet at this angle as an orientation 1 
current sheet. At t = 47.8 mins, we see that a new 
current sheet has formed at a similar location, but 
now at an angle of approximately 37r/4 relative to 
the positive ^-direction. We shall refer to a cur- 
rent sheet formed at this angle as an orientation 2 
current sheet. At t = 49.2 mins, we see that a cur- 
rent sheet has formed again in a similar location, 
but that this current sheet is again of orientation 
1. Finally, at t = 50.8 mins, we see the formation 
of an orientation 2 current sheet, again in a simi- 
lar location (it is also clear that this current sheet 
is weaker, i.e. |J| m ax is decreased, compared to 
the previous current sheets). Thus, Figure 3 illus- 
trates the formation of a cycle of current sheets, 
i.e. the formation of orientation 1, followed by 
orientation 2, followed by orientation 1 again, and 
so on. Note that this figure is a qualitative illus- 
tration of the periodic nature of the current sheet 
formation in this system (§3.2 will provide quan- 



titative evidence). Note that our terminology of 
orientation 1 and 2 is purely arbitrary; a simi- 
lar periodic cycle of current formation was seen in 
McLaughlin et al. ( 2009) and was referred to a 
cycle of horizontal and vertical current sheets. 

Thus, we recover the results of Murray et al. ( 
2009). Murray et al. demonstrated that, using 
fieldline tracing, one can quantitatively demon- 
strate the periodic change in connectivity of the 
open flux over time, i.e. evidence of reconnection. 
We recover the fieldline-tracing results of Murray 
et al. (their Fig. 4). Hence, we use the ter- 
minology oscillatory reconnection to refer to the 
periodic formation of orientation 1 and 2 current 
sheets as well as the associated periodic changes 
in connectivity. We now focus on the observable 
consequences and outputs of such an evolving sys- 
tem. 

3.1. Generation of quasi-periodic outflows 

As current sheets form and reconnection com- 
mences in the simulation, we observe strong out- 
flows emanating from the ends of the current 
sheet, in agreement with classic steady-state re- 
connection theory (e.g. Sweet 1958; Parker 1957; 
Petschek 1964). Upon leaving the current sheet, 
these jets (strong outflows) collide with the mag- 
netic field already in the outflow regions and are 
deflected into two secondary jets at angles of ap- 
proximately ±7r/4 to the original jet (a termina- 
tion shock is also present). The schematic struc- 
ture of these reconnection jets is in good agree- 
ment with the description of Forbes ( 1988). Since 
these secondary jets are deflected at angles of ap- 
proximately ±7r/4 and due to the orientation of 
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Fig. 4. — Contours of (top row) v x and (bottom row) v y , in units of km/s, and selection of fieldlines at times 
t = 46.7, 47.8, 49.2 and 50.8 mins. Blue/red corresponds to positive/negative motions. White dot indicates 
the fixed point (xo, yo) — (— 2.5Mm, 7Mm). Dotted line denotes a change from equilibrium transition region 
to coronal temperature profile. 



the current sheets, they manifest themselves as ei- 
ther primarily horizontal and/or vertical outflow 
jets. For orientation 1 current sheets, jets from 
the lower end of the current sheet give rise (pe- 
riodically) to negatively-directed v x and v y mo- 
tion. Meanwhile, jets from the upper end of the 
current sheet give rise (periodically) to positively- 
directed v x and v y motion. For orientation 2 cur- 
rent sheets, the two jets have mixed velocity com- 
ponents, and this can be clearly seen in Figure 
4. Figure 4 shows contours of v x (top row) and 
v y (bottom row) at four snapshots in our simu- 
lation (the blue/red colour table corresponds to 
positive/negative motions) . 

Let us first consider the v x behaviour (top row 
of Figure 4). The characteristic behaviour can be 
summarised as follows: 

• At t = 46.7 mins (orientation 1 current sheet), 
we have strong horizontal outflows, with positively- 
directed v x motions (blue) ejected from the up- 
per end of the current sheet. 

• At t = 47.8 mins (orientation 2 current sheet), 



we have negatively-directed v x motions ejected 
from the upper end of the current sheet. 

• This cycle repeats, and we have positively- 
directed v x motions again at t = 49.2 mins 
(orientation f) followed by negatively-directed 
v x motions at t = 50.8 mins (orientation 2). 

A similar pattern is observed emanating from the 
lower end of the current sheet but with the op- 
posite orientation, i.e. at t = 46.7 mins (ori- 
entation 1) we have negatively-directed v x mo- 
tions, followed by positively-directed v x motions 
at t = 47.8 mins (orientation 2), again in a repeat- 
ing cycle. 

Let us now consider the v y behaviour (bottom 
row of Figure 4). The characteristic behaviour can 
be summarised as follows: 

• At t = 46.7 mins (orientation 1 current sheet), 
we have strong vertical outflows, with positively- 
directed v y motions (blue) ejected from the 
upper end of the current sheet. Interestingly, 
we also have negatively-directed (but weaker) 
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Fig. 5. — (a) Time evolution of v x (km/s) at the point (xo, yo) = (—2.5 Mm, 7 Mm). (6) Time evolution of 
v y (km/s) at (xo,yo) = (—2.5 Mm, 7 Mm). The red dashed lines indicate an exponentially damped envelope 
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where A = 0.1464s . In both figures, the dotted line denotes zero velocity. 



v y motion adjacent to (to the left of) our 
positively-directed v y outflow. This negatively- 
directed v y motion is associated with the inflow 
region of the simulated current sheet, and thus 
acts at right-angles to the current sheet orien- 
tation. 

• At t = 47.8 mins (orientation 2), we now have 
a positively-directed v y motions ejected from 
the upper end of the current sheet. Again, 
negatively-directed v y motion appears adjacent 
to (now to the right of) our positively-directed 
v y outflow (again associated with inflow into our 
reconnection region). 

• The cycle then repeats: we have positively- 
directed v y motions again at t — 49.2 mins (ori- 
entation 1) and negatively-directed v x motions 
at t — 50.8 mins (orientation 2). 

The vertical v y motions are of obvious interest for 
comparison with observations. Following the ter- 
minology of Murray et al. ( 2009) we refer to these 
vertical v y motions as a collimated jet. 

3.1.1. Transverse/ swaying collimated jets 

We now consider these vertical outflow jets and 
swaying, transverse motions in further detail. To 
do so, we measure the v x and v y signal at a fixed 
point (xo,yo) — (—2.5 Mm, 7 Mm) where this par- 
ticular point has been chosen as it is located close 
to the upper end of the orientation 1 current sheet 
(slightly above and to the right) and thus is well 



placed to measure the transverse motions and v y 
outflows of the resultant jets. 

By considering the evolution of the colli- 
mated jet, we note that the central axis will 
be horizontally-displaced periodically as the cur- 
rent sheet contracts and lengthens (i.e. evolves 
between orientation 1 and 2, and back again). 
Whilst evolving from orientation 1 to 2, the col- 
limated jet will appear to move in the negative 
x— direction (i.e. an orientation 1 current sheet 
first contracts, then lengthens into an orientation 
2 current sheet). Conversely, whilst evolving from 
an orientation 2 current sheet (back) to orien- 
tation 1, the collimated jet will appear to move 
in the positive x— direction. This displacement 
will repeat itself as the cycle repeats. Thus, the 
collimated jet displays a characteristic swaying 
or transverse motion. Note that this transverse 
behaviour is specifically due to the oscillatory re- 
connection mechanism, and would be absent for a 
single, steady-state reconnection jet. 

Qualitatively, this transverse displacement can 
be seen by comparing the locations of the (blue) 
collimated jet in Figure 4 (bottom row). How- 
ever, such displacement can also be quantitatively 
measured from our simulation. In Figure 5a, we 
see the evolution of v x (km/s) at the fixed point 
(xo,Vo) = (-2.5 Mm, 7 Mm). The oscillatory be- 
haviour of v x can be clearly seen, i.e. the quasi- 
periodic transverse/swaying displacement. Note 
that the change from orientation 1 to 2 occurs at 
t w 47 mins. Before this time, the strong positive 



v x motion is associated with the initial formation 
of the orientation 1 current sheet. 

3.1.2. Quasi-periodic vertical outflows 

Figure 5b shows the evolution of v y , i.e. the 
vertical outflow, at the fixed point (x ,j/o) = 
(—2.5 Mm, 7 Mm). We can clearly see the oscilla- 
tory behaviour: the outflow changes from v y < 
at t s=s 45.5 mins, to v y > at t w 46.5 mins, to 
v y < at t w 47.5 mins, to v y < at t w 49 mins. 
After £ rs 50 mins, the signal is much weaker since 
the bulk of the (oscillatory) reconnection has oc- 
curred by this time. 

We fit the oscillation with an exponentially- 
damped envelope ~ Vj,-rr 



-At 



-At 



for v y > and 
for v y < (red dashed lines), where 
A = 0.1464 s _1 is determined experimentally. It 
is also important to note that there is a stronger 
positive v y signal than negative v y signal, i.e. 
there is more 'upflow' than 'downflow'. This can 
be most clearly seen by calculating the average 
of the two exponentially decaying envelopes, i.e. 
~ 0.5 (vy-max + %_min) e~ At , and this average is 
denoted by the blue dashed line. Since this curve 
is positive, it indicates the presence of a positive 
upflow, i.e. for zero flow, we would expect the 
two exponentially decaying envelopes to average 
to zero. 

Note that the generation of this positively- 
directed vertical net flow is tightly localised to 
the upper ends of the (evolving) current sheet, 
and thus one must be careful when analysing the 
oscillatory signal at a fixed point (as in Figure 
5b). Similar oscillatory results and upflows are 
obtained for other choices of fixed points located 
around (#o>2/o) = (—2.5 Mm, 7 Mm) and within 
the same (north-easterly) domain of connectivity. 

We also note that for both Figures 5a and 5b, 
the conditions before t w 43 mins are explained by 
the initial flux emergence itself, i.e. the flux tube 
is expanding upwards and sideways into the solar 
atmosphere and so, for xq < 0, a small v y > and 
v x < is initially present. 

3.2. Quantitative analysis of periodicity 

In §3.1, we reported on the formation of strong, 
quasi-periodic, vertical outflows generated by os- 
cillatory reconnection. As detailed in McLaughlin 
et al. ( 2009), by considering the evolution of the 




Fig. 6. — Time evolution of A z located at the cen- 
ter of our (moving) current sheet. Insert shows 
a blow-up over a shorter time period. Vertical 
dotted lines denote how we measure the period of 
oscillation. 



vector potential at the center of the current sheet, 
one can obtain a quantitative measure of the os- 
cillatory nature of the system. In our numerical 
experiment, we track the center of the (rising) cur- 
rent sheet and plot the value of the vector poten- 
tial at that point versus time. This evolution can 
be seen Figure 6. 

For 45 < t < 60 mins, the oscillatory nature 
of our system is clear, and changes in the vec- 
tor potential can be related to changes in connec- 
tivity of the system (see §4 of McLaughlin et al. 
2009). Meanwhile, for 60 < t < 70 mins, the vec- 
tor potential changes very little and then increases 
(70 < t < 90 mins) to a plateau (t > 90 mins) to 
a constant A z = 5.527. The oscillations between 
45 < t < 60 mins are directly related to the oscilla- 
tory reconnection mechanism at work in our simu- 
lations (i.e. a local effect, searching for an equilib- 
rium) whereas the behaviour for 60 < t < 70 mins 
and beyond are related to the system settling into 
equilibrium on a global scale. There is no change 
in the vector potential for t > 90 mins. 

Since we are primarily interested in oscillatory 
reconnection in this paper (as opposed to the gen- 
eral behaviour of flux emergence) we shall focus on 
the time period 45 < t < 60 mins, and a blow-up 
of the evolution of A z over this time is shown in 
the insert of Figure 6. We measure the period of 
oscillation by considering two local extrema in the 
vector potential (indicated by the vertical dashed 
lines in the insert). We measure a period of 195 
seconds, i.e. 3.25 mins, for the parameters and 
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Fig. 7. — (a) Parametric study of the evolution of A z located at the center of the current sheet for eight 
values ranging from 2 < Bq < 3, i.e. 2.6 x 10 3 < B < 3.9 x 10 3 G. Note that the colour bar indicates the 
value of the non-dimensional Bq . (b) Parametric study of the evolution of A z for Bq — 1.9 (failed emergence, 
green line), Bq = 3.5 (blue line), Bq = 4 (purple) and Bq = 5 (yellow). Bq = 3 is plotted as a dashed line 
to aid comparison between both subfigures. 



initial conditions we have considered. 

4. Parameter Study of Bq 

We now investigate the dependence of our 
system to varying Bq, i.e. the initial magnetic 
strength of our buoyant flux tube. This will allow 
us to see how the period of oscillation changes and 
investigate how robust our system is to the onset 
of oscillatory reconnection. 

Figure 7a shows the full time evolution of A z 
(located at the center of the current sheet) for sev- 
eral different values of the initial value of the mag- 
netic field strength of the buoyant flux tube, i.e. 
we vary the value of Bq. Since this is a numeri- 
cal parameter study, we choose to present values 
associated with Bq, i.e. the non-dimensionalised 
initial magnetic field strength. Specifically, Fig- 
ure 7a reports the evolution of eight values rang- 
ing from 2.6 x 10 3 < B < 3.9 x 10 3 G or, in 
non-dimensional units, 2 < Bq < 3. Note that 
§3 above corresponds to Bq = 3.25 x 10 3 G, i.e. 
Bq = 2.5. From the resultant curves, it is clear 
that oscillatory behaviour is present in all of these 
numerical experiments. Thus, we can measure the 
corresponding period of oscillation for each choice 
of Bq and this is shown in Figure 8a. 

We also investigate larger and smaller values 
of Bq and these results can be seen in Figure 7b, 



which presents values of Bq = 2.47, 4.55, 5.2 and 
6.5 x 10 3 G, i.e. Bq* = 1.9, 3.5, 4 and 5. These 
choices yield behaviors of a significantly different 
nature to those seen in Figure 7a. This is for 
two reasons: firstly, for values of Bq < 2, i.e. 
Bq < 2.6 x 10 3 G, the initially-submerged flux tube 
does not successfully emerge into the solar atmo- 
sphere. This is in agreement with the results of 
Murray & Hood ( 2006) who found that for low 
initial magnetic field strengths, the tube cannot 
fully emerge into the atmosphere since the buoy- 
ancy instability criterion is not satisfied. Thus, in 
this model, for values of Bq < 2.6 x 10 3 G, we have 
no flux emergence (or failed emergence) and thus 
no onset of oscillatory reconnection. Thus, high 
in the solar atmosphere, A z remains constant. In 
Figure 7b, we plot A z corresponding to a Bq = 1.9 
failed emergence (constant value green line) for 
comparison. 

Figure 7b also shows the time evolution of the 
vector potentials corresponding to Bq = 3.5 (blue 
line), Bq = 4 (purple) and Bq = 5 (yellow). We 
note that the evolution of A z is again significantly 
different to that seen in Figure 7a. Bq = 3 is plot- 
ted as a dashed line to aid comparison between the 
subfigures. The reason for this significant change 
in behaviour is that for Bq > 3, we have successful 
flux emergence and current sheet/ X-point forma- 
tion in a similar manner to that detailed in S3 



10 



(b) 

location of X-poirt at time=100 mins 
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Fig. 8. — (a) Parametric study of period of oscillation for eight values ranging from 2 < Bq < 3, i.e. non- 
dimensional Bq. (b) Parametric study of the location of the X-point after t = 100 mins for eight values 
ranging from 2 < Bq < 3. Note that the colour bar indicates the value of the non-dimensional Bq and the 
choice of axes are only a subset of the total numerical domain. 



above but, critically, these strong current sheets 
now eject plasmoids from their ends and this fun- 
damentally changes the properties of the current 
sheet /X-point, since the plasmoids take magnetic 
flux with them as they leave the current sheet /X- 
point. Thus, although we still have oscillatory be- 
haviour present in the behaviour of A z correspond- 
ing to Bq =3.5 and 4, this represents a fundamen- 
tally different regime to that seen for 2 < Bq < 3 
in Figure 7a. The effect is even more pronounced 
for Bq = 5 and beyond. 

Let us now investigate the period associated 
with our choice of Bq (Figure 8a) where we re- 
strict ourselves to 2 < Bq < 3 (as explained above, 
values above and below these limits correspond to 
fundamentally different regimes). We find that the 
period of oscillation increases with Bq , for periods 
105 — 212.5 seconds, i.e. 1.75 — 3.5 min. 

Finally, we investigate how the final location of 
the X-point depends upon our choice of Bq. We 
find that the location X-point at t = 100 mins, i.e. 
the final/resting location, is displaced both higher 
and further to the left for stronger values of Bq 
(again we restrict ourselves to 2 < Bq < 3). This 
behaviour can be seen in Figure 8b. 

5. Discussion and Conclusions 

We have performed numerical experiments of 
magnetic flux emerging into a coronal hole, mod- 



eled as a pre-existing unipolar magnetic field, 
within a stratified solar atmosphere, and solve 
the compressible and resistive MHD equations us- 
ing a Lagrangian remap, shock capturing code: 
LARE2D. The long-term evolution of the system 
is followed and we investigate how the reconnect- 
ing magnetic systems behave as they search for 
an equilibrium. We find that the initial rise and 
expansion of the emerging magnetic flux tube are 
in good agreement with that reported in the exist- 
ing literature (Figure 2). Reconnection is initiated 
at t = 41.6 mins (for Bq = 2.5 system) in which 
inflows acting perpendicular to the current sheet 
bring magnetic field into reconnection region, and 
magnetic flux is ejected at the ends of the current 
sheet. We find that oscillatory reconnection occurs 
in the model (Figure 3) i.e. a process in which re- 
sistive diffusion at the X-point is coupled to global 
advection of the outer field. We find that the first 
reconnection reversal, i.e. change from an orien- 
tation 1 current sheet to an orientation 2, occurs 
around t w 47 mins. The mechanism for oscilla- 
tory reconnection is well described by McLaugh- 
lin et al. ( 2009) and Murray et al. ( 2009) and 
occurs due to a local imbalance of forces, primar- 
ily the gradients in thermal pressure, between the 
neighbouring flux systems. 

We find strong horizontal outflows from both 
ends of current sheet and, once oscillatory recon- 
nection is initiated, these change direction peri- 
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odically. Similarly, we find strong vertical out- 
flows in the positive y— direction emanating from 
the upper end of our current sheet, and these 
sit side-by-side with negative y— direction inflows 
bringing magnetic flux into the reconnection re- 
gion. The direction of these positive/negative ve- 
locities changes as the system changes from orien- 
tation f to 2 (Figure 4) but we also find on average 
there is a vertical upflow directed in the positive 
y— direction (Figure 5b). 

We find that the vertical outflows/collimated 
jet will be displaced horizontal during the shorten- 
ing and lengthening of the evolving currents sheets 
which results from the oscillatory reconnection 
mechanism. This displacement gives the verti- 
cal outflow jets a swaying nature and explains the 
transverse nature of the oscillations. This was fur- 
ther confirmed by analysing n I (-2.5Mm, 7Mm, t), 
which clearly demonstrated the transverse be- 
haviour (Figure 5a). 

In order to quantitatively estimate the vertical 
outflows, we measured the vertical velocity at a 
fixed point: v y (— 2. 5 Mm, 7Mm, t). Here, the pos- 
itive (material ejected from the end of the current 
sheet) and negative (inflow into the reconnection 
region) vertical flows were visible and as well as 
the periodic nature of the oscillations (Figure 5b) . 
Only a few complete periods were observed and 
so we labelled this quasi-periodic behaviour. We 
fit the damped oscillations with exponential en- 
velopes and found that there was a preference for 
positive velocity over negative, i.e. evidence of a 
positively-directed vertical net flow in the model. 
Note that the generation of this positively-directed 
vertical net flow is tightly localised to the upper 
ends of the (evolving) current sheet and so simi- 
lar results are obtained only for other choices of 
fixed points located around (—2.5 Mm, 7 Mm) and 
within the same domain of connectivity. 

By tracking the vector potential, A z , at the 
center of the (moving) current sheet, we quantita- 
tively measure the period of oscillation, allowing 
us to measure a period of 195 seconds (3.25 mins) 
for an initial magnetic flux tube strength of B = 
3.25 x f0 3 G (B% = 2.5). 

We also perform a parameter study varying the 
initial magnetic strength of the buoyant flux tube, 
i.e. Bq. We find that for a range of parameters ap- 
plicable to the solar atmosphere, 2.6 x 10 3 < Bq < 
3.9 x 10 3 G, we observed similar flux emergence 



and oscillatory reconnection behaviour, with each 
Bo corresponding to its own period of oscillation 
in the range 1.75 < period < 3.5 mins. Essen- 
tially, the stronger the initial flux tube strength, 
the longer the period of oscillation. 

However, for B Q < 2.47 x 10 3 G, we do not 
observe successful flux emergence and thus, ob- 
viously, there is no subsequent oscillatory mo- 
tion. This is in agreement with Murray & Hood 
( 2006) who found that for low initial magnetic 
field strengths, the tube cannot fully emerge into 
the atmosphere since the buoyancy instability cri- 
terion is not satisfied, i.e. 'failed' flux emergence. 
Further details of the buoyancy instability crite- 
rion can be found in Newcomb ( 1961), Yu ( 1965), 
Thomas & Nye ( 1975), Acheson ( 1979) and Ar- 
contis et al. ( 2004). Thus, it is the buoyancy in- 
stability criterion that dictates the lower limit in 
our model, given our choices of parameters. For 
different parameters, e.g. a stronger/weaker equi- 
librium uni-directional magnetic field, the buoy- 
ancy instability criterion will have a higher/lower 
threshold, although a full investigation must be 
undertaken to investigate the true behaviour. 

We also investigate larger values, i.e. Bq > 
4.55 x 10 3 G. Here, unlike for 2.6 x 10 3 < B Q < 
3.9 x 10 3 G, we observe the formation of plas- 
moids ejecting from the ends of our current sheet. 
These ejected plasmoids change the properties of 
the X-point, e.g. taking magnetic flux with them. 
Thus, even though we still have oscillatory be- 
haviour, seen in the evolution of A z , this repre- 
sents a fundamentally different regime than that 
of 2.6 x 10 3 <B < 3.9 x 10 3 G. The exact reason 
for plasmoid formation above a particular thresh- 
old strength is uncertain, and will be investigated 
in future work. 

Finally, we investigated how the final location 
of the X-point depends upon our choice of B n . We 
found that the location X-point at t = 100 mins, 
i.e. the final/resting location, was displaced both 
higher and further to the left for stronger values of 
Bq (we restricted our records to 2.6 x 10 3 < B < 
3.9 x 10 3 G) . Both the longer duration periods of os- 
cillation and increased displacement/height of the 
X-point can be fully explained since stronger B 
tubes have larger emergence velocities and thus 
greater momenta. Thus, the larger momentum 
flux tubes are able to compress the X-point to a 
greater extent (resulting in a stronger/longer cur- 
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rent sheer and thus longer periods for the subse- 
quent oscillatory reconnection, which explains Fig- 
ure 8a) and, secondly, higher Bo flux tubes with 
larger momenta carry the tube higher and farther 
into the atmosphere (explaining Figure 8b). 

Thus, we have presented numerical simulations 
that naturally generate quasi-periodic flows with 
a characteristic transverse/swaying aspect. Such 
outputs result from the oscillatory reconnection 
physical mechanism, i.e. in a self-consistent man- 
ner since no periodic driver is imposed on our sys- 
tem. The vertical speeds of the outflows, 20 — 
60 km/s, arc comparable to those reported in re- 
cent observations. By varying the initial strength 
of our submerged flux tube, we recover periodic- 
ities in the range of 1.75 — 3.5 mins. Thus, the 
mechanism of oscillatory reconnection may pro- 
vide a physical explanation for the generation of 
some of the recent quasi-periodic, transverse, ver- 
tical motions/jets/outflows reported by a multi- 
tude of authors (see §1 for details). In particular, 
the oscillatory reconnection mechanism presented 
here may explain the observations and simulations 
by Nishizuka ct al. ( 2008; 2011) and He et al. 
( 2009a; 2009b), who detail Hmode/SOT obser- 
vations and simulations of transverse motions on 
a spicule, originating from the cusp of an inverted 
Y-shaped structure, as well as more recent work 
by Ding et al. ( 2011) and Harra ct al. ( 2011). 

The oscillatory mechanism presented here may 
also partially explain quasi-periodic pulsations 
(see reviews by Aschwanden 2003; Nakariakov 
& Melnikov 2009). Such oscillatory behavior has 
been observed in a number of solar and stellar 
flares (Mathioudakis et al. 2003; 2006; McAteer 
et al. 2005; Inglis et al. 2008; 2009; Nakari- 
akov et al. 2010; 2011) but the true physical 
mechanism responsible remains uncertain. 

Finally, it should also be noted that the mech- 
anism generates both v x and v y motions together 
(one cannot exist without the other), and such mo- 
tions are exponentially damped. Thus, if such sig- 
nals are detected, they may be decaying not due to 
a damping mechanism but due to the generation 
mechanism itself . This is not surprising given that 
flux emergence injects a finite amount of energy 
into the (oscillatory reconnection) mechanism, and 
so it is expected that the phenomena and outputs 
will also only be of a finite duration, i.e. this is a 
dynamic reconnection phenomenon as opposed to 



the classical steady-state (time-independent) re- 
connection models. 

Although we have only presented a specific ex- 
ample of reconnection initiated by flux emergence, 
we believe the oscillatory reconnection mecha- 
nism described in this paper is a robust, general 
phenomenon that may be observed in other sys- 
tems that demonstrate finite-duration reconnec- 
tion. Further studies should focus on the inclusion 
of heat conduction (e.g. Miyagoshi & Yokoyama 
2003; 2004) which is expected to reduce the tem- 
perature of the outflow jets. However, the density 
of outflows will also be increased by heat conduc- 
tion, i.e. to ensure force balance in the current 
sheet, which may make the outflow jet more ob- 
servable (e.g. Shiota et al. 2005). Finally, the os- 
cillatory reconnection mechanism itself should be 
investigated further (e.g. Gruszecki et al. 2011) 
and extended to fully three-dimensional studies. 
Evidence of oscillatory reconnection in 3D flux 
emergence simulations was recently reported by 
Archontis et al. ( 2010). 
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